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ABSTRACT 


The author employs a 16K linked, FORTRAN II version of the 
problem oriented PACER simulation program written by Prof. Paul 
T. Shannon, Dartmouth College, Hanover, New Hampshire to 
develop a non-linked, FORTRAN IV program which is compatible with 


the GE 625 system installed at the University of Kansas. 


The author then formulates a steam power system model and 
tests several configurations of this system to demonstrate 

(1) the efficacy of the modified PACER procram and 

(2) the potential usefulness of such a modified PACER 

program in either graduate or undergraduate chemical 


engineering curriculum. 
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Chapter I 


INTRODUCTION 


General. In early 1965 the National Science Foundation estab- 
lished a project at the University of Michigan to evaluate the current 
and future uses of computers in engineering design. education. The 
final report issued by this project observed that the majority of 
analytical work done in upper-level engineering courses was ofa 
routine nature, where the solution procedures are generally brief and 
straightforward and where there is normally only one correct answer. 
For this reason, the report stated, engineering design courses had 
fallen into general disregard end in some instances they have been 
eliminated from the curricula. Rather than take a pessimistic view from 
these circumstances, the report concluded that: 

Such courses can, with the introduction of the computer, 

take on a new luster and make significant contributions 

to the education of our engineering graduates, In the 

past, the engineering student has only had time to 

investigate perhaps one or two possible (and usually 

not very effective) designs. With the new hardware and 

software systems, the student designer will for the first 


time have the opportunity to investigate a large number 
of feasible solutions to his design problem.... 


In effect, because of his exposure to a large number of 
designs and his ability to manipulate system configurations 
and design variables easily, the student can gain some 
facets of "engineering experience" which cannot now be 
acquired outside the industrial environment (and which 
would require a much longer elapsed time). Thus, in 
addition to its use as an analysis tool, the time shared 
computer with appropriate design-oriented software 
should enable the student to improve his intuition and 
understanding of the nature of good design: this is the 
raison d'etre of engineering design courses Qi 


One of the most promising approaches to the demand for efíeciive 
engineering design education lies in the recent development of problem- 
oriented approaches to programming languages. The problem-oriented 
approach invoives the development of a general processor or executive 
program which is capable of handling a general class of problems with 


a minimal amount of programming effort on the part of the problem 





designer. This frees the designer to concentrate on the primary 


objective of system design and analysis. 


Of the problem-oriented simulation programs developed thus Íar, 
one of the most notable is PACER. PACER was developed by Prof. Paul 
T. Shannon of Dartmouth College and was developed primarily to 
process material and energy balances of typical system design problems 
encountered in chemical engineering. Solutions are processed by a 
modular or unit operaticns concept with iterative calculations under- 


taken when recycle conditions are encountered. 


Objective and Scope of Thesis. This thesis will examine the 
historical development of digital simulation programs culminating with 
the development and some applications of PACER, The essential 
concepts and characteristics of PACER as developed by Prof. Shannon 
will be described and illustrated. The author will further discuss the 
modifications made to PACER to make it compatible with the computation 
system installed at the University of Kansas. Additionally, the author 
will trace the development of a semi-realistic model of a steam power 
plant developed to evaluate the modified PACER program. As a related 
objective, the usefulness of the model in an undergraduate chemical 


engineering education curriculum will be explored. 


Finally, the advantages and limitations of the modified PACER 
program and suggestions for its further improvement and refinement 


will be presented for consideration. 





Chapter II 
HISTORICAL DEVELOPMENT OF DIGITAL 
SIMULATION PROGRAMS AND 


APPLICATIONS TO ENGINEERING EDUCATION 


Historical Development of Digital Simulation Programs. Initial 
attempts toward computer simulation of chemical engineering systems, 


were made in the early 1950's and were limited almost solely to analog 
computers. The digital computer at that time required tedious machine- 
language programming and it was slow and expensive for such uses. 
Accordingly, it could not effectively compete with the analog Somme 
(Figure 1 presents a schematic illustration of the historical development 


of simulation programs). 


R. G. Selfridge is generally credited with the development, іп 
1955, of the first digital simulation program for use on the IBM 701 
comente d This initial effort was followed during the period from 1956 
to 1958 by the development of a group of digital computers known as 
digital differential analyzers (D.D.A.'s). The D.D.A,'s were designed 
to duplicate the functioning of the analog computer, but proved to be 


exceedingly slow and EE 


One of the most significant advances in digital simulation 
occurred with the development of procedure- oriented languages, such 
as FORTRAN, in the late 1950's. The availability of these languages 
facilitated the heretofore cumbersome prograrnming task by permitting 
the programmer to write his program in terms closely resembling the 
symbolic representation of proklem variables and parameters. A com- 
piler converted the symbolic program into machine language instructions. 
The FORTRAN language additionally was provided with a library of 
functions and subroutines which performed commonly encountered mathe- 
matical operations such as logarithims, exponentation, etc. FORTRAN 
coupled with increasingly spectacular improvements in memory capa- 
bility and speed of operations, where operations are currently specified 
in nano-seconds (billionths of a second), gradually begen to assert the 


effectiveness of digital simulation programs. 
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Figure 1. Development of Simulation Programs? 
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The reiative advantages and disadvantages of digital and analog 
simulation became more distinct as each system was progressively 
developed. It was readily acknowledged that the analog computer 
enjoyed the advantage of parallel operation which obviated the careful 
attention to logical execution sequence and extensive anticipation of 
alternatives required by tne sequential nature of the digital computer. 
Digital computers normally require extensive reprogramming for different 
problem structures. Mathematical integration, a natural operation for 
an analog computer, requires specific and extensive programming on a 
digital computer. On the other hand, the digital computer provides 
certain specific advantages which the analog computer cannot approach: 
greater reliability, memory and logic capabilities as well as a decidedly 


better dynamic range. 6 


An extensive effort was undertaken to incorporate the best 
features of the analog and digital computer into a digital simulation 
capability referred to as a digital-analog simulator. These simulation 
systems generally consist of block-oriented’ subroutines which perform 
or simulate the natural functions of an analog computer such as 


integration and differentation. 


Two representative examples of digital-analog simulators are the 
MIDAS and MIMIC programs. The most significant capability of MIDAS/ 
is its integration function which automatically adjusts the integration 
step to a specified size and combines this feature with a fifth-order 
Milne, predict-correct integration routine. These characteristics 
provide the user with assurance that the integration is being performed 


correctly. 8 


The next effort in digital simulation involved the development of 
equation-oriented, digital-analog simulation programs which did not 
rely upon the block-oriented structure imposed on MIDAS and its 
contemporaries. Chief among this group of programs were DS1/90? 
and MIMIC. 19 


MIMIC, introduced in 1965 by the authors of MIDAS, differs in 
concept from the less flexible block-oriented MIDAS. MIMIC offers 


several other features not comparably available in MIDAS, suchas 
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a wider variety of available subrovtines, FORTRAN comoatibility and 


intuitive solution capability for implicit algebraic equations. ы 


A separate class o£ programming languages has been designed to 
simulate systems characterized by a sequence of discrete events such 
as schedule optimization and resources allocation. Representative 
examples of this type of program are GASP, 12 GES, 13 SIMSCRIPT, 13 
еы: 


Тһе most significant current advance in the field of digital 
simulation is the development of problem-oriented or application- 
oriented programming languages which are designed to handle a specific 
class or specific classes of problems. Some of the more significant of 
this type of program are: 

14 15 


1. BLODIB ` and ECAP ` for signal analysis and circuit 
design in electrical engineering 
Zu: pyana? © for vibration analysis in mechanical engineering, 


and 
Mm" 19 | | | 
= PACER for system simulation in chemical 


engineering. 


PACER, an acronym for Process Assembly Case Evaluator Routine, 
was conceived in 1961 by Prof. Paul T. Shannon while he was in the 
employ of Humble Oil and Refining Company. After extensive modifi- 
cation and development, PACER received its first practical test in 
October 1964 when Prof. Shannon joined with Prof. A, I. Johnson of 
McMaster University, Hamilton, Ontario, Canada and Dr. Norman 
Cooke from Canada Industries Ltd. to undertake a simulation of C-I-L's 
300 ton/day contact sulfuric acid plant in Hamilton. This initial under- 
taking was successfully completed and reported in April 1965. The 
sulfuric acid plant simulation involved the use of 42 equipment, infor- 
mation gatherinc or parameter setting subroutines which were inter- 
connected by approximately 70 information transmittal streams. Тһе 
model thus conceived involves the solution of approximately 500 
simultaneous equations involving almost 1000 stream variables and 200 
equipment parameters. A complete simulation run required 8 minutes of 


computer time on the IBM 7040 computer. 





In 1965 the staff and the senior chemical engineering students 
at McMaster University joined with C-I-L in the design of a 600 ton/day 
sulfuric acid plant based on metallurgical off-gas. The results of this 
study were presented to C-I-L and representatives of other companies 
in April 1966. Today, more than twelve industrial firms are using 
PACER in the analysis of their processing problems, and thirty-nine 
universities are integrating PACER into their teaching and research 


activities. 


Use of Computers and Simulation Programs in Engineering 
Education. In 1959 the College of Engineering of the University of 
Michigan undertook, under the sponsorship of the Ford Foundation, a 
study of the feasibility of broad scale integration of electronic computer 
use into the educational process. The final report of this study 
recommended that: 

(1) All engineering students be introduced to computing 

techniques early in their college education, preferably 

at the freshman or sophomore level. 

(2) That they develop a proficiency in a procedure- 


oriented language such as FORTRAN, ALGOL, MAD, etc., 
and 


(3) Required upper-level engineering courses incorporate 

the requirement for computer solutions of several engineering 

problems of graded difficulty by the student prior to gradu- 

ation. 

A disciplinary appendix to the final report entitled "The Use of 
Computers in Chemical Engineering Education" cited the purely 
educational aspects of computers in the education process, stressing 
the salutary effects of precise formulation, definition, and logical 


organization required for computer solution of engineering problems. 


A survey of engineering school deans and department heads 
taken in 1962 as a part of the Project revealed that approximately 30% 
of all engineering students attended schools which had a required 
introductory computer course for undergraduate students. The study 
further indicated 45% of engineering students attended colleges or 
universities which required students in some engineering disciplines 
to attend such an introductory course. A subsequent study conducted 
during the 1965-56 academic year showed that approximately 87% of 


all engineering undergraduate students received some required 





instruction in digital computer techniques, An analysis of this data 
indicates that in the near future almost all undergraduate engineering 
students will attend a required digital computer introductory course. *! 
Quality, however, apparently has not matched quantity as the 
summary report on "Computers in Engineering Design Education" 
concluded. (See Chapter I). Most of the computer solutions assigned 
in upper level engineering courses are routine in nature, One major 
shortcoming in the field of engineering design has been traditional 
reticence of design instructors to employ sophisticated mathematical 
analysis techniques in their courses. Further, existing computer 
hardware and programming languages, due primerily to their general 


nature, were not amenable to flexible use in design problems. 


However, two recent developments in the field of computer 
technology hold the promise of a vast improvement in the area of 
engineering design education. The first of these new developments is 
the time-shared computer with remote terminal installations. These 
facilities will permit a large number of users to have almost instan- 
taneous access to the computer and will permit these users to communi- 
cate with the computer in an almost conversational manner. This time- 
sharing capability, coupled with the development of a new hierarchy of 
problem-oriented computer languages and processors such as PACER, 
permits the visualization of the investigation of a large number of 
feasible configurations of complex engineering problems by a student 
in the course of a semester. In effect, the computer will then take 


f 22 
on the role of an effective imparter of "engineering experience". 
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Chapter III 


CONCEPT, ORGANIZATION AND OPERATIONS OF PACER 


Concept of PACER, As previously mentioned PACER was developed 
by Prof. Shannon to provide a highly-flexible, systematic method for 
the simulation and design of complex chemical engineering systems and 
processes. The basic concepts underlying the PACER system may be 


summarized as follows: 


1. PACER is a systematic method which permits a progressive 
development of comprehensive models of chemical engineering systems 


or processes. 


2. The processing or operational sequence and the physical 
configuration oí the system or process are included in the input data 
provided by the user and can be readily modified to explore alternate 


system or parametric configurations. 


3. Mathematical models of individual equipmenis are available 


for use during development of the total system or process. 


4. Earlier work and results may be incorporated into the process 


model with minimal effort on the part of the user. 


5. The successful completion of each analysis contributes to 
the development of a librarv of mathematical equipment models which 


à ; 23 
may be used in subsequent analyses. 


PACER conceptually then differs from a procedure-oriented 
language such as FORTRAN, FORTRAN permits a problem solver to 
write solutions for mathematical and engineering problems in a language 
which symbolically approximates commonly-encountered engineering 
and mathematical terms. It does not provide a general or even a specific 
solution method. It's greatest advantage lies in its flexibility--the 
ability to represent many types of problems, ranging from the social 


sciences and law to nuclear physics, with the same language. 


This generality also represents a deficiency of FORTRAN. Its 
use, in many cases, requires the complete restructuring of complex 


solution programs for only minor system or probiem configurations, 
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Such a requirement can be exceedingly cumbersome to a systems 
designer or analyst who is primarily interested in the results of the 


simulation. 


PACER was deveioped to provide to a chemical engineer the 
ability to rapidly and accurately solve varying configurations of complex 
chemical engineering systems. PACER accomplishes this by providing 
a rigidly-structured program which is capable, with minimal effort on 
the part of the user, of solving complex chemical engineering systems. 
No tedious rewriting of the program is necessary for changes in config- 
uration. PACER, then, provides the logical next step in computer 


language design-- Ñmacro-procedure-oriented programs. 


It permits a user to systematically describe (by the use of a 
highly-structured data deck) the topology and parametric constraints 
of the system to be solved, Major changes in the actual physical 
configuration oí a system can be accommodated with relatively simple 
changes to the data deck. It also provides the user the option of 
specifying or not specifying the order in which equipment will be cal- 
culated. If the order of calculation is not specified, PACER enters into 


a trial and error procedure to determine the solution order. 


Ihe PACER concept of operations may be effectively illustrated 
by the presentation of a simple problem involving (1) a flow mixer, 
(2) a heater and (3) a flow splitter. Figure 2 presents the flow diagram 
for such a system, as well as the major components of the data deck 
which would be used to communicate the system configuration and 
parametric representation to the PACER executive routine. A more 
extensive discussion of the data deck and its components is presented 


subsequently in this chapter and Appendix A. 


Organization of PACER, The major elements of PACER are: 


l. A mainline program which serves as the general processor 


for PACER operations and conducts al! iterative calculations. 


2. Information modifying and parameter setting subroutines. 
These include: 
(a) EQCAL - Calls proper equipment subroutines to calculate 


system or process equipments. 
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b. Process Matrix 
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с. Equipment Parameters 
STREAM 
NG) = | STREAM VARIABLES 
l; 0 ES, 80. TO 0 100. 0 
D. 0 15; 60. 200. 200. 0 0 
Зи 0 lo; 0 0 0 0 0 
4, 0 15. 0 0 0 0 0 
m 0 oe Q 0 0 0 Ü 
6. 0 5, 0 0 0 0 0 


d. Stream Variables 
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e. Convergence Criteria 
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(b) TEST - if calculation is in an iterative (recycle) mode this 
subroutine tests the values of the output stream variables for the equip- 
ment being calculated to determine if the calculated values meet con- 


vergence criteria established by the user. 


(c) XDATA - provides the space for additional equipment 


variables if required by tne user. 


(d) MAJOR - called by PACER when order of calculation of 
equipment is not supplied by user. MAJOR directs search for list of 
equipments which will calculate as output streams up to three assumed 
known input streams. When such à list is developed the equipments are 
calculated as either direct or iterative calculations depending on their 


contiguration. 


(e) EQUIP - serves as a bookkeeping subroutine for all equip- 
ments which are calculated directly, Bookkeeping functions for equipment 


calculated in an iterative mode are performed by the mainline program. 


(f) GUESS1 - forms under the direction of MAJOR an ordered 
list of unknown input streams. From this list, up to three unknown 
streams are assumed known and a search is made to find a list of equip- 
ment which will calculate as output streams the assumed known input 


streams. 


(4) GUESS2 - carries out under the direction of GUESS1 the 
search for a list of equipment which will calculate as output streams 
those unknown streams assumed known in GUESS1. Ií such a list of 


streams can be found, GUESS2 calls for appropriate calculation. 


(h) DREAD - reads data deck for problem configuration (run). 


A new set of data is required for each configuration (run). 


(i) EXREAD - reads additional equipment parameters if pro- 


vided by user as well as order of calculation if specified. 


0) SEFLAG - sets numerical indicators (flags) indicating 
status of equipments and streams provided as data. Flags may indicate 
either that an equipment is used or not used or that a stream is 4 feed 


stream, product stream, internal stream or not used. 


(k) DPRINT - prints data supplied by user. Primarily used as 


a debugging feature. 
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(1) EXPRNT - prints additional equipment parameters supplied 


or calculated by PACER as well as order of calculation specified by user. 


(m) RPRINT - prints results of run conducted by PACER. Results 


are printed at completion of each run. 


3. Equipment subroutines which are mathematical models of 
equipments to be calculated. PACER as supplied by Prof. Shannon con- 
tained forty-two equipment subroutines in the PACER equipment subrou- 
tine library. The user may supply such additional equipment subroutines 


as he may require. 


4. Data deck - the data deck describes the system or process 
configuration and any known stream variables or equipment parameters. 
It also includes certain dimension and control parameters as well as 
certain other operational indicators or controllers. The principal con- 
stituents of the data deck are: 

(a) Process Matrix - the process matrix is the principal 
segment of the data deck. It describes the sytem or problem configuration 
and indicates stream interrelationships. Column 1 contains the equip- 
ment subroutine which performs the mathematical calculation of the equip- 
ment listed in column 1. The columns after column 2 indicate the numbers 
of the streams flowing into and out of equipment number listed in column 
1. All input streams are listed first. All output streams must be preceded 


ужа нише Sign. 


(b) Equipment Parameters Matrix - contains for each equipment 
Jisted in the process matrix a list of equipment parameters such as weight, 
volume, power, or heat transfer coefficient. Column 1 indicates the 
equipment number. All columns after column 1 contain equipment para- 
meters for that equipment. The equipment parameters are characteristic 


of that equipinent only. 


(c) Stream Variables Matrix - contains stream variables for 
each stream listed in the process matrix. The stream varíables represent 
sucn characteristic items as pressure, temperature, specific heat, totel 
flow rate, component flow rates, etc. Stream variables are consistent 
for each stream. Column 1 contains the stream number. Column 2 is a 


flag indicating the status of the stream. If the flag equals 1, the stream 
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is a feed stream. If it equals 2, it is a product stream. If it is equal 
to 0, itis an internal stream. The remaining columns contain the values 


oí tne stream variables for that stream. 


(d) Convergence Criteria - contains in each column an abso- 
Jute fractional convergence value for each corresponding column in the 
stream matrix. Convergence of an item in the stream variables matrix 
occurs when the value calculated for that variable does not vary from the 
current value of that variable stored in the stream variables matrix by 
more than the absolute convergence criterion. The convergence criterion 


for a particular variable location (column) is the same for all streams. 


(e) Calculation Order Listing - if provided or specified by the 
user it indicates the order in which the equipment will be calculated as 
well as the mode (direct or iterative) by which the equipments will be 
calculated. Column l indicates the equipment number and column 2 
the mode of calculation. A value of 1 in the calculation mode column 
indicates that the equipment is to be calculated directly. A value of 2 
indicates that the equipment is to be calculated iteratively. If the cal- 
culation mode is 3 the equipment is the last equipment in a list of 


equipment to be calculated in an iterative mode. 


(f) Equipment Subroutine Listing - a listing which establishes 
the correspondence between the subroutine name and subroutine number 


to facilitate calling of proper subroutine for equipment calculations. 


(g) Other control and informational parameters and data such 
as problem dimension parameters, problem control parameters, and title 


of simulation run. These items are specifically definea in Appendix A. 
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Chapter IV 


MODIFICATIONS MADE TO PACER 


Modifications Made to PACER, The PACER source deck supplied 
by Proí. Shannon was impressed on a 2400 foot magnetic tape at 556 
bits/inch, 80 characters/record, blocked 1. With the assistance of 
the University of Kansas Computation Center, the magnetic tape was 
converted into an interpreted, punched card deck of approximately 
Унион оке: 


Detailed analysis of the interpreted deck revealed that: 


(1) The program was organized into eight 16K links as well as 
a starter link and several sample problems. The eight links were 
arranged as follows: 

a. LINK 1 - contained mainline segment for conducting iterative 
calculations, calling and testing subroutines and certain equipment 
subroutines used in the contact sulfuric acid plant simulation. 

b. LINK 3 - consisted of mainline seginent which directs search 
for list of equipment when order of calculation is not specified, sub- 
routines for formulation of equipment calculation list, and a subroutine 
for recording values calculated in direct calculations. 

C. LINK 4 - composed of a mainline segment which directed 
input/output operations and subroutines which effected input/output 
Operations and stream and equipment flagging. 

d. LINKS 2,5,6,7,8 - contained remaining equipment and 
information modification and gathering subroutines used in contact 
sulioric acid plant. 

(2) The linking procedure employed in the program was not 
compatible with the GE 625 system installed at the University of 
Kansas. 

(3) The program was written in FORTRAN II, 


Further analysis indicated that the 11,500 cards could be 
reduced to approximately 2800 by removing all those cards pertaining 
to equipment subroutines used in the contact sulfuric-acid-plant 
simulation as well as those used in the saniple problems at the end of 


the source deck. The cards remaining after this modification constituted 
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the PACER executive routine. However, it stil] existed in a non- 
compatible linked version and was written in FORTRAN II which was 


not permitted on the GE 625 system. 


Accordingly, two additional modifications were required to 
make the PACER executive routine compatible with the GE 625 system. 
First, the program had to be converted to a linked version acceptable 
by the GE 625 and the program had to be rewritten in FORTRAN IV, A 
detailed investigation of the organization of the executive routine and 
the simulation model planned for development indicated that the PACER 
executive and the simulation model could be combined into a non- linked 


program which would not exceed the memory available on the 625. 


The conversion from a linked to a non-linked version created 
several exacting problems. First, the combination of the linked main- 
line segments into one non-linked mainline program created a statement 
number redundancy. Second, it caused the invalidation of a major 
portion of the linked program transfer logic. This problem required 


extensive effort to resolve. 


The rewriting of the executive routine in FORTRAN IV required 
considerably less effort than the conversion to a non-linked program. 
The major changes required were the rewriting of arithmetic function 


statements and input/output instructions. 
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Chapter V 


DEVELOPMENT OF STEAM SYSTEM MODEL AND 


ТЕСИ ТЕ ОБИЛИЕ PACER EXECUTIVE ROUTINE 


Objectives of Testing. The principal objectives of testing 
conducted with the modified PACER executive routine were: 


(1) To verify the validity of the modifications made to the 
original version of PACER; 

(2) To demonstrate the capability of PACER to simply and 
accurately simulate the performance of varying configurations of com- 
plex chemical engineering systems. In this regard it should be 
pointed out that specific care was taken throughout the course of the 
modification of PACER to insure that no alterations were made in 
either PACER logic or method of solution. Finally, 

(3) To provide the experimental and factual basis for eveluating 
the potential usefulness of the modified PACER executive routine in a 
chemical engineering curriculum. This objective was facilitated by 
the steam power plant model developed by the author for testing with 


the modified PACER program. 


Development of Simulation Model, Feasibility or validity testing 
of the modified PACER program presented the contrasting alternatives 
of (1) simulation of a segment of the contact sulfuric acid plant by 
use of the equipment subroutines supplied by Prof. Shannon or (2) 
development of a new simulation model. The second alternative was 
selected since it provided the opportunity for development of a more 


thorough understanding of PACER program logic. 


The simulation model selected was that of a conventional 
steam power system utilizing fuel oil as its energy source. The reasons 


supporting this selection were: 


(1) The familiarity of the author with steam power systems. 
This familiarity was gained by the author from service for two years 
as an instructor in Thermodynemics at the United States Naval 
Academy, as well as practical acquaintance with steam power systems 


as the primary means for naval ship propulsion. 
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(2) The general applicability of steam systems to chemical 
engineering. 
(3) The ready adaptability of such a model to the undergraduate 


chemical engineering curriculum at the University of Kansas. 


Accordingly, the author undertook the development of the 
essential equipment subroutines necessary for the approximation of 
a Steam power system. As envisioned, the steam power system would 


consist of: 


(№) 2s Turbine 

(2) A Condenser 

(3) A Cooling Water Circulating Pump 
(4) A Condensate Purnp 


(5 A Heater Mixer 
) A Feed Pump 
) An Economizer 
) A Boiler 
) A Flow Splitter 
(10) A Superheater 
) A Fuel Oil Heater 
) A Forced Draft Blower 
) A Flow Mixer 


Such a system, it was envisioned, would be amenable to demon- 
strating the ability of PACER to perform iterative material and energy 
balances by the multiple application of individual equipment sub- 


routines. 


The basic system as envisioned and developed by the author is 
presented in Figure 3. Figures 4 and 5 present alternate configurations 
of the basic system. Table 1 provides a stream legend for Figures 3-5. 
Figure 6 presents the process diagrams for the basic system on 
enthalpy--entropy coordinates. Analysis and development of mathe- 
matical models for individual components of the basic system are 


presented in Figures 7-19. 
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Table 1. 


да 
Description 
Steam 
Exhaust Steam 
Impurities 
Losses 
Condensate, 
Cooling Water 
Cooling Water 
Cooling Water 
Auxiliary Steam 
Auxiliary Steam 
Condensate 
Auxiliary Steam 
Auxiliary Steam 
Feed 
Feed 
Auxiliary Steam 
Auxiliary Steam 
Feed 
Steam 
Fuel Oil 
Boiler Fuel Oil 
Boiler Fuel Oil 


Boiler Blowdown of Impurities and 
Accompanying Feed 


Fuel Oil 
Superheater Fuel Oil 
Superheater Fuel Oil 
Auxiliary Steam 
Auxiliary Steam 

Air 

Air 

Auxiliary Steam 


Auxiliary Steam 


Steam Power Plant Stream Legend 
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Stream Number Description 
33 Flue Gases 
34 Steain 
35 | шск Б 
36 Make-Up Feed 


Table 1l. (continued) 
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Power = MN x (h, - հշ) x No 
12 ü (hy И һу) е m 

s un - steam flow rate 

հլ - entering enthalpy 

No - internal engine efficiency 

Ма - mechanical efficiency 

Үй - isentropic expansion enthalpy 

h, - exhaust enthalpy 


Figure 7. Turbine 





Q steam * (h, 
W ow - ЗИ Ж (r; 
W. =F, x W 
inpur 1 steam 
loss ` E, Է W steam 
= W вн + W - W 
cond steam inpur muf loss 
CHELAN = ТЕЕ = ШК ЕТА? 
LN իր Ը 
foe BAZ 
ARFA = ОЛСО x BHETAM) 
ET "рап У ess 
О - heat transferred 
411 - inpurities flow rate 
W - make-up feed flow rate 
muf 
T; Ts - temperature 
Իշ - losses flow factor 
ТА - water flow rate in blowdown 
lu. - system losses 
THETAM - log mean temperature difference 
THETAl - larger temperature difference 
THETA2 - smaller temperature difference 
U - overall heat transfer area 
A - condenser heat transfer area 
w - steam flow rate 
steam 
h „Ва, ће, де - enthalpy 
Ce - cooling water specific heat 
Py - inpurities flow factor 
W ‚ — condensate flow rate 
сопа 


Figure 8. Condenser 
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DN ON 
Ա. - saturated steam flow rate 
Wa - pumped liquid flow rate 
H - Specific volume of pumped fluid 
Pe: Po ~ pressure 
h,,hig - enthalpy 
hos - isentropic expansion enthalpy 
NSi - internal engine efficiency 
Nr - mechanical efficiency 


Figure 9. Circulating Pump 
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(hy 9 Ë hy3) x Ni, 


fl 5) 


= Np (hyo хо 
= he + V x (P,1 7 Ps) 
W - saturated steam flow 
aux 
V - specific weight of pumped fluid 
T - pressure 


һул, уз - enthalpy 

ћ) - isentropic expansion enthalpy 

N. - internal engine efficiency 

М - mechanical efficiency 

Мя - weight rate flow of fluid being pumped 


Figure 10. Condensate Pump 
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- feed flow rate 
Мз! Ура Мод: Уза 


- total energy into heater-mixer 


- saturated steam flow rates 
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energy supplied to condensate 


h - enthalpy 
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11 temperature 


- specific heat of condensate 


- condensate flow rate 


Figure 11. Неаїег- Міхег 
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5717 
Шы x Bent) 
(hy g hi) x N. 
Ly = бе“ Веб“ бб Х Ма 
п б = h), + V x (Pis = Ру) 
W - saturated steam flow 
aux 
Wa - pumped liquid flow rate 
V - Specific volume of pumped fluid 
Pis: P14 - pressure 
һуб, 17, Вуд - enthalpy 
ћ 65 - isentropic expansion enthalpy 
Nr - internal engine efficiency 
N. - mechanical efficiency 


Figure 12. Feed Pump 
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Q - heat transferred 
Wig - feed flow rate 
C £a - specific heat of feed 


Түр 115 - feed temperature 
T 


145 - flue gas temperature 


33: 
ра - flue gas flow rate 
= - flue gas specific heat 


THETAM - log mean temperature difference 


THETAl - larger temperature difference 


THETA2 - smaller temperature difference 
U - overall heat transfer coefficient 
AREA - economizer heat transfer area 


Figure 13, Economizer 
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Wyo] = QO/(HV x N.) 
: ТҮЙЕ 
НУ - heating value of fuel oil 


C,H,O,S - carbon, hydrogen, oxygen and sulfur weight 
fractions in fuel oil 


uen - total boiler blowdown 

a E inpurities flow rate 
Pa - inpurities blowdown factor 
pr - pure water in blowdown 
M - saturated steam : 
Q - heat transferred 
ћ 1 - enthalpy of saturated liquid 
Па, По - enthalpy 
W -fuel oil flow rate 
N, - boiler efficiency 


Figure 14. Boiler 
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No machemaical calculations take place in 
the splitter. The splitter merely serves as 

a convenient manner for splitting the sat- 
urated steam flow produced in the boiler into 
its various components without needlessly en- 


cumbering the boiler calculations. 


Figure 15. DS NUNC 1 
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steam is (h, y hg 4) 


14600 ИО 62000 x (H - 0,/8) ООо 


= Q/(HV x №) 
О - heat transferred 
հլ Р h34 - enthaipy 
HV - heating value of fuel Ji 
С,Н,0,,5 - carbon, hydrogen, oxygen and sulfur weight 
fractions in fuel oil 
W - fuel oil flow rate 
Ғо2 
М, - boiler efficiency 


Figure 16. Superheater 
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Q - heat transferred 
Ire – fuel oil flow rate to boiler 
Wo - fuel oil flow rate to superheater 
C £o - Specific heat of fuel oil 
t5 4 Ct, g) - fuel oil exit temperature for proper atomic dica 
11 ( էշ 6) - fuel oil entrance temperature 
— saturated steam flow rate 
aux 
Пот, Под - enthalpy 


Figure 17. Fuel Oil Heater 
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A/F - air-fuel ratio 
XA - percent excess air 


Ф 


stoichiometric coefficient for oxygen 


1 
Co - stoichiometric coefficient for nitrogen 
С; - stoichiometric coefficient for fuel oil 
JO) - molecular weight of fuel oil 
А - air flow rate 
ա о? - fuel oil flow rates 
W - saturated steam flow rate 
aux 
V - velocity 
g - acceleration due to gravity 
J - Joules constant 
հշլ,հշշ - enthalpy 


ey 
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isentropic expansion enthalpy 


31s 
N. - mechanical efficiency 
Nos - internal engine efficiency 


Figure 18. Forced Draft Blower 
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5 - flue gases flow rate 
М, Мс - fuel oil flow rates 
Wap - air flow rate 


Figure 19. Mixer 
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Plan of Testing. The development of a series of approximations.to 
achieve the objectives established above required a systematic analysis 
to determine if any of the objectives were complementary or interdepend- 


ent, and, if so, in what manner and in what order. 


An analysis of these relationships revealed that the objective of 
verifying the validity of the modifications made to PACER was the objec- 
tive on which all the other objectives depended. Additionally, it was 
concluded that objectives two and three could be achieved as derivatives 
of objective one. That is, if simulation tests could be developed which 
would establish the validity of the modifications and the accuracy of the 
model, then the efficacy of PACER. and the model and their usefulness in 
chemical engineering education could be determined by derivation from 


those results. 


Accordingly, a series of six simulations were scheduled for testing 
with the modified PACER program, The first, Run 01, was designed to 
approximate the performance of the basic system with the order of calcu- 
lations specified by the user. The second simulation, designated Run 03, 
involved the same basic system but did not specify the order of calcula- 
tions. The accuracy of the solutions was not of primary concern in these 
tests. Rather, the objective was to ascertain if the modified program were 
processing the data in accordance with the methods described in the docu- 
mentation provided by Prof. Shannon and inferred from an analysis of the 
Original PATER, 


Two additional system configurations were tested to demonstrate the 
ability of PACER to provide accurate solutions of varying configurations of 
the same system. The first of the altered configurations, labeled Run 02, 
involved all the components of the basic system except the superheater. 
The second altered configuration consisted of the basic system less the 
economizer. This run was designated Run 04. Both Runs 02 and 04 were 
run with the order of calculation specified and unspecified. The results 
obtained from these runs were compared with hand-calculated values and 
corrections as necessary, and were made in the models until the computed 


values agreed with the hand calcuations for all cases. 


The final results of all runs are presented in Appendix B, with 


extensive diagnostic print-out (a PACER option) provided for Runs 01 and 





ՀՄ 


03. A guide for the interpretation of the results of the simulations accom- 


panies the results in Appendix B. 
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Chapter VI 


ANALYSIS OF TEST RESULTS 


Accuracy of Results. Acceptable limits of accuracy vary with the 
nature of the probiem being evaluated, Quite clearly, the acceptable 
limits of accuracy in determining the configuration of the critical mass 
of a fissionable material would be much more rigorous than those to be 
imposed on the sytem parameters and configuration of a steam power 
plant, In the latter case, slide rule accuracy is normally acceptable. 
Accordingly, a convergence criterion of .02 (296)--an approximation of 
slide rule accuracy--was selected for all steam variables other than those 
non-varying items such as stream number and flag value. With this cri- 


terion, allruns converged to solution in no more than five iterations. 


The values of stream variable calculated by PACER were compared 
to hand-calculated values and in all cases were within slide rule accuracy 


of hand-calculated values. 


Effect of Specifying Order of Calculations. PACER undertakes trial 





and error (iterative) solutions in all cases where the order of calculations 
is not specified or where an iterative calculation is specified. Under 
these conditons, the advantage of specifying the order of equipment cal- 
culations was not clear. In Run 02 (without superheating) both the calcu- 
lations with and without the order of calculations specified converged in 
five loops. In Run 04 (without economizer) the calculetions with the order 
of calculations specified converged in three iterative loops, whereas the 
calculations without the order of calculations specified converged in four 
loops. Runs 01 and 03 (basic system) converged in five and four loops 


respectively. 


The inconclusive nature of these results would seem to indicate 
that the PACER method of determination of order of calculation is as efíec- 
tive in promoting convergence as intuitive designation of the order of 


calculation--or at least no more ineffective. 


Another aspect of this consideration is the comparison of the values 
of stream variables obtained with different orders of calculations. Quite 


clearly convergence values will be expecied to vary depending on the 
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logical location o£ equipment in the order of calculations. While this 
expectation was confirmed, differences were in all cases less than the 


convergence criterion. 


Effects of Changes in System Configuration. Two major changes 
to the basic steam power system were simulated with the modified PACER 
program, The first change involved removal of the superheater from the 
basic system. This change necessitated the use of saturated steam in 
the turbine with the result that for the same turbine power requirement, 
steam flow increased 12%. However, fuel oil flow increased only 1%. 


Condenser surface area increased 6%. 


The secona change involved removal of the economizer from the 
steam power system. This change necessitated the addition of the heat 
that would have been added by flue gases in the economizer as the heat 
of combustion of additional fuel oil. This change resulted in a 20% 


ШШерсасе ІП fuel oil flow. 


Realism of the Models. The qualitative and quantitative nature of 
the stream variables tended in all cases to verify the realistic nature of 
the models developed by the author. For example, the use of saturated 
steam in lieu of superheated steam would undoubtedly require more steam 
to produce the same turbine power. The greater steam flow would norm- 
ally require more fuel oil. However, the fact that the additional steam 
was saturated and the original superheating load was no longer required 


had an offsetting effect on fuel oil flow. 


Additionally, as was demonstrated in Run 04, the loss of the regen- 
erative heating affect of the economizer would quite clearly require more 


fuel oil for the same power requirement. 


While these results tend to verify the realistic nature of the steam 
power system developed by the author, it should be understood that no 
attempt was made in the development of the model to completely depict 
reality. Rather an attempt was made to provide a semi-realistic model 
which could be refined subsequently. It should be noted in this regard 
that most of the data and mathematical bases for equipment and sytems 
models are based on reliable data developed from many years in the 


desion, development and operation of naval steam power plants. 
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Validity of Modifications Made to PACER, The close conformance 
of the values obtained in the simulations with hand-calculated values, 
as well as the general and specific correspondence of results to format 
and arrangement of data as described in documentation provided by Prof. 
Shannon clearly establishes the validity of modifications to the original 


PACER program. 





Chapter VII 


AN ASSESSMENT OF PACER 


General. Tne author has demonstrated the concept, principles and 
method of application of the PACER simulation program by the develop- 
ment of a complex model of a steam power system, and successfully 
approximating its performance on a modified version of the original PACER. 
The validity of the modifications made to PACER have been conclusively 
established by the accuracy of the results achieved. The realistic nature 
of the model has been tentatively established by the qualitative and 


quantitative nature of the results. 


PACER As A Simulation Device. PACER has been demonstrated to be 
a powerful, flexible and efficient method for simulation of complex con- 


figurations of chemical engineering systems. 


However, certain areas where improvements in PACER method of 
operations or concept could be made were noted and are discussed here 


for edification of potential users. 


(1) The addition or removal of equipment subroutines from PACER 
is cumbersome with present arrangement of calling segments and sub- 
routines and LPHA listing. A more permanent, non-mnenonic name for 
equipment subroutines requiring only a change in subroutine name to con- 
form to permanently established but unused name is one method of allevi- 
ating this problem. This would obviate changes in mainline program which 


become more difficult when executive program is on tape. 


(2) Consideration should be given to establishment of permanent 
STRMI and STRMO matrices to accommodate all streams on a permanent 
basis for a run, This would facilitate multiple calling of stream variable 
and information transmittal in a particular run. However, it could develop 
that the ultimate result would be to trade the present inflexibility of a 


temporary STRMI and STRMO matrix for an inflexible permanent matrix. 


(3) Consideration should be made for incorporating capital invest- 
ment costs into equipment subroutines. This information could take the 
form of a cost which varies with heat transfer area, pressure, temperature, 


etc. Incorporation of this data into the basic structure of an equipment 


= 





me 
| yh 

subroutine is a relatively simple matter-- primarily a function of availa- 

bility of cost information and its relationship to design or operating 


parameters. 


- Use of PACER In A Chemical Engineering Curriculum. Іп evaluating 


ЕЕ О спре simulation of the C- I- L contact sulfuric acid plant 

using the PACER program, Prof. Shannon remarked: 
It should be emphasized that this study, although it 
involved a unique and very close cooperation with 
industry, was in the main an exercise in the education 
of engineering students, The students were concerned 
with the detailed modelling of specific equipment, yet 
had a very broad exposure to a real industrial plant and 
its over-all behavior. They learned of the difficulties 
and importance of obtaining good plant data, and of 
writing successful programs of greater complexity than 
previously encountered, The project was admittedly a 
very ambitious one; however the large effort expended 


has placed PACER into the hands of the staff for future 
use in the undergraduate and graduate program. 


The staff generally agree that most of the students have 

benefitted by this exposure to the quite new concept of 27 

information handling and decision making utilizing PACER, 

PACER then has already demonstrated its effectiveness in instilling 
in students some aspects cf "engineering experience" which would not 
normally be obtainable outside an industrial environment, This, as 
stated in the University of Michigan Report on the uses of computers in 


engineering education, is the raison d'etre of engineering design courses. 


It is further considered that PACER can be an effective means of 
developing a systematic approach to problem solving by chemical engi- 
neering students, The demands of computer programming in development 
of equipment subroutines and systems concepts requires strict attention 


to problem structuring and formulation. 


Recommmendations. In consideration of the results of the simu- 
lations conducted with the modified PACER program and the conclusions 


drawn above, it is recommended that: 


(1) Consideration be given to the incorporation of the moditied 
PACER program into existing undergraduate courses involving basic prob- 


lems in chemical engineering. The steam power system as developed by 


28 
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the author is amenable to immediate use on an individual or term project 


basis, and 


(2) After suitable evaluation of the effectiveness of the initial use 
of PACER that consideration be given to further and more extensive use of 


PACER in graduate and research projects at the University of Kensas. 





EO; 


11. 


12, 
13. 


14. 
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Аррепдіх А 


USER'S GUIDE FOR THE MODIFIED PACER SIMULATION PROGRAM 


General. The most significant feature of the PACER simulation 
method is its ability to perform realistic and accurate simulations of 
varying configurations of complex chemical engineering systems with a 
minimal effort on the part of the user. Yhe element of PACER which 
makes this capability possible is a highly-structured data deck which 
permits the user to fully characterize system configuration and the re- 


straints or parameters under which it operates, 


Successful use of PACER is'dependent upon strict adherence to a 
prescribed set of instructions for the formulation, assembly and modifi- 
cation of the data deck, It is the purpose of this appendix to provide 
the necessary instructions in a manner which will permit a user to 
employ PACER's powerful calculational skilis without having to master 
either the logic or structure of PACER. One exception exists to this 
condition, That occurs when it is necessary to develop or modify equip- 
ment subroutines or add them to or delete them from the program. In that 
case a minimal knowledge of the structure and logic of PACER is required 
to insure the subroutines conform to PACER logic and that they are properly 
entered into PACER structure. Both of these aspects will be treated in 


this appendix, 


System Conceptualization. The initial step in preparing a system 
configuration for simulation by PACER is to establish the system concept 
and commit the concept to formalization by preparation of a flow diagram. 
This flow diagram should break equipment assemblies into subassemblies 
and equipments capable of being mathematically modelled. For example, 
a boiler normally contains an economizer and a superheater as integral 
equipments, and when one normally refers to a boiler it is assumed that 
he includes the economizer and superheater if they are provided. However, 
breaking these elements out and handling them as separate equipments 
facilitates mathematical and system modelling as well as increasing the 
flexibility of the model. Further, the functional components if of suffi- 
cient general nature could be used in multiple applications within a 


single simulation or as a related equipment in subsequent simulations. 
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In developing the system flow diagram all equipments, the per- 
formance of which will be simulated by mathematical models (equipment 
subroutines), should be numbered. No redundancy in equipment number- 
ing is permitted although equipment numbers need not be successive. 
While equipments are being numbered, they should also by named. For 
convenience it would be preferable to assign them a six-character alpha- 
numeric name which is mnemonically related to the equipment function, 
e. g. TURBNS for turbine and CNDNRS for condenser, No redundancy is 
permitted in the last 3 characters of a unit name, although individual 
characters within a unit name may be repeated. For example, TURBNS 
and TURBN6 are acceptable, while SPLTRS and SPHTRS are not. 


After all equipments have been numbered and named, all inter- 
connecting streams should be numbered. No redundancy is permitted 
in stream numbers. However, stream numbers need not be successive 


nor is redundancy between stream and equipment numbers prohibited. 


Eguipment Subroutine Format. The equipment subroutine is the one 
facet of PACER that poses a programming burden on the user. If an 
equipment subroutine which satisfies the analytical and calculational 
requirements of an equipment which is to be simulated is not available 
in the PACER subroutine library, the user must develop an equipment 


subroutine which accomplishes these objectives. 


Several challenges confront the user developing an equipment sub- 
routine for use in a PACER simulation. The first challenge is to insure 
that the equipment to be modelled is, in fact, the equipment which should 
be modelled. In this respect a general effort should be made to break 
equipment assemblies down into their functional elements, that is, as 
previously mentioned, a boiler into a boiler, an economizer and a super- 
heater. Mathematical analysis is facilitated and simplified by this 
approach, Further, specific effort should be directed toward incorporating 
generality into subroutines. That is, an economizer is in essence a two- 
fluid heat exchanger. If the equipment subroutine can be written fora 
general two-fluid heat exchanger which is in this application called an 
economizer, then the equipment subroutine can be used in subsequent 
future applications where a two-fluid heat exchanger is employed. The 
fact of this generality should be documented into the subroutine by 


appropriate comments. 





ЭЭ 
Once the mathematical representation of an equipment has been 
developed it is necessary to translate this into the syntax employed by 
PACER. This task can be mastered quickly by knowledge of a few basic 
concepts of PACER. 


PACER performs equipment calculations by using mathematical 
relationships provided by the user in the equipment subroutine and 
stream and equipment data provided by the user in the stream (SN) and 
equipment (EN) matrices of the data deck (to be discussed further in 
section dealing with the data deck). PACER calls data out of the stream 
matrix into temporary STRMI (inlet streams) and STRMO (outlet streams) 
for the equipment to be calculated. The inlet and outlet streams for a 
. particular equipment are designated by the user in the process matrix 
(KPM). Stream variable information is then ready for use in PACER cal- 
culations, Equipment parameters are called from and returned directly 
to the EN matrix for the equipment to be calculated. Once the equipment 
has been calculated the values in the STRMI and STRMO matrices are 
returned to permanent storage in the SN matrix. The STRMI and STRMO 
matrices are then available for temporary storage of inlet and outlet 


stream variables for next equipment to be calculated. 


It is then incumbent upon the user to transcribe his mathematical 
notation into a notation which indicates to PACER where in the STRMI, 
STRMO or EN matrices the variable or parameter to be used is to be lo- 
cated. For example, let us convert the mathematical equation for deter- 
mining turbine stream flow rate into a notation usable by PACER. The 


equation (neglecting dimensional constants) is: 
W = P — Wk 
In FORTRAN this is represented as 
W = P/WK 
However, in PACER notation it would be written as 


STRMI (1,10) = EN (NE,3) / (STRMI(1,5) - STRMO(1,5)) 
(where WK = (STRMI(1,5) - STRMO(1,5))) 


In the above notation 
STRMI(1,10) - indicates the tenth variable in the first inlet stream 


to the equipment being calculated. It also represents the tenth 





56 
variable in the stream matrix vector for the stream which is the 
first input stream to the equipment. In the context of this sample, 


this variable represents the stream flow for that stream. 


EN (NE,3) - indicates the third parameter for equipment NE. Equip- 
ment NE is the equipment number for the equipment being calculated 
and will be set elsewhere in PACER. The turbine power rating for 


equipment NE will be stored in this location. 


STRMI (1,5) and STRMO (1,5) - represent the fifth variable in the 
first inlet and outlet streams respectively. In this example these 
variables take on the value of the entering and leaving enthalpies 


of the first inlet and outlet streams to equipment NE, 


The effect of the above calculation is to calculate a value for 
steam flow for the stream number which is the first inlet stream to equip- 
ment NE, After calculation this calculated value is stored in the SN 
matrix releasing the STRMI and STRMO matrices available for further 
calculations. Since the values of EN(NE,3) and STRMI(1,5) and STRMO 
(1,5) were not modified, the corresponding values in the EN and SN 


matrices will not be modified. 


Figures Al. and A2. are the equipment subroutines for the mathe- 
matical models developed for the turbine and condenser in Chapter 4 (See 
Figures 7 and 8). It should be pointed out that the assignment of variable 
and parameter locations in the SN, EN, STRMI and STRMO matrices are 
the responsibility of the user. Further, stream variable location must 
be consistent for all streams in the SN, STRMI and STRMO matrices. 
That is location 5 in the SN, STRMI and STRMO matrices will throughout 
the course of a run represent the same variable, but not necessarily the 
same value. For example, in the illustrative example discussed above, 
location 5 represented the enthalpy for all streams but the enthalpy for 
all streams was not necessarily the same value and some streams did not 
even store a value for enthalpy if it were not relevant to the calculation 
of that equipment. In such a case PACER normally assigns a value of 


zero to that location. 


One additional note of caution is that a common and dimension 


deck identical to that used in the mainline program must be the first item 
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. Co*oe aPACER COHMON AND DIMENSION DFC: 
D]I"ENSION STRMI(B,15), STRHCI (5,9), STRMO(B8,15), STRNCO(5,5), 


24. 


SUPROUTINE TURRNS 


1S5N(50,15)4, SNC(50,5), ENCS04 19) 1 ENC C50, 5) , МНЕСА | (53), МАМЕ (50), 
2NA"E2(30) , NAMPES(30Y , NAHEA (50) , KEt (C3U ) , KSEM(S0, 3) , KUHSC BD y , KPSCÓ), 


ЗОЕ! SCA5), DELSC(5, *PMC3SO0,9) ,KEFLAGEGO),KSFLAGGCOQ), TITLED 4), 2 
&LPHAC30) ,KOMPNCO((10),AEN(C4, 49) 4 NELISTOCSO,2) , RSPRNT(C16),PAPCR(ES,19), 
SNPAPER( 20) ' ' 
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COMMON KPM, NYMAX, NIMAX, TITLE ° ВЕ ЕН 


COrMON AEN, LPHA, KOMPNO, NOCCMP, NELIST 
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COYMON LIMIT, LIMIT2, LOOP, KRACK, KBACK2, KVOID 
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COrMON MCHAIN, NCHAIN, KLOCKi 

CO-MON KOUNT1, KOUNY?, KOUNTS, KOUNT4; KOUNTIS 37557-27: 


CO“MON EN, МЕ, NLIST, LIST 


CO" MON  XUS,NUSMAX, KPS, NPS ELE : me - 
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Сечевеб |. НСИТ | МЕ СМОМЧ5 DETERMINES THE COOLING WATER FLOW RATE AND THE SyQ- 
FATE AREA REQUIRED TO CONDENSE THE EXHAUST STEAM FROM THE TURBINE; 

IT ALSO CALCULATES THE FLOW PATE OF INPURITIES INTRODUCED INTO THE 
SYSTEM AND THE STEAM AUD WATER LASSES ву LEAKAGE, EVAPORATION,- ETC, 
TH= COWDENSER SURFACE AREA [S CarCULATED FPOM THE RELATIONSHIP THAT 

С z U X A X THETAM, “HERE 11 TS THE OVERALL HEAT TRANSFER COEFFICIENT 
AN: THETAM IS THE LOG MEAN TEMPERATURE DIFFERENCE. 
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NO"ENMCLATURE 
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Figure A2. Subroutine CEONDES 
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in each subroutine. This facilitates communication of information from 


subroutines to mainline and other subroutines. 


Incorporation or Removal of Equipment Subroutines into PACER 
Structure. The tlexible nature of PACER permits the incorporation or 


removal of equipment subroutines to conform to core memory limitations, 
to permit the simulation of completely different systems or to permit 
modification of existing systems. Such flexibility is not achieved, how- 
ever, without. a commensurate demand on the user. In this case, it is 
the requirement that the incorporation or removal of a subroutine requires 


a change in the semi-permanent structure of PACER, 


Each subroutine, in addition to being assigned a six-character, 
alpha-numeric title, is assigned an equipment subroutine number depending 
on the location of its calling statement in PACER. Equipment subroutines 
1-10 are called from EQCAL and subroutines 11-30 are called from the 
mainline program. The location of these calling statements and maximum 
number of subroutines are a function of the system to be simulated. How- 
ever, a note of caution. An increase or decrease in the maximum number 
of equipment subroutines or any other dimension establishing parameter 
requires a change in the appropriate matrix or vector in the common and 


dimension deck. 
The incorporation of a subroutine requires 


(1) Obtaining an equipment subroutine number from these not 


currently being used. 
(2) Writing a calling statement for that subroutine. 
(3) Providing appropriate transfer directions or transfer indicators. 


For example, the addition of an equipment subroutine COOLR8 as 
subroutine number 5 would require the insertion of the following FORTRAN 
cards in EQCAL immediately following statement number 4 and its transfer 
card, if any, and immediately preceding statement number 6: 

Column Noro 7 
9 CALL COOLRS 
RETURN 
The insertion of the same subroutine as equipment subroutine number 15 


would require the insertion of the following cards in the mainline program 
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immediately following statement number 14 and its following transfer 


card: 
Column No. 4 7 = 
15 CALL COOLR8 
GORO 77/7 


Or the addition of COOLK8 as equipment subroutine number 21 or any 
number greater than 21 would require the insertion of the following cards 
in the proper location following statement number 5201: 
Column No. 4 7 
ZEN TIAS OO0LRO 
GO TO 7777 
The removal of ЩЕ equipment subroutine requires the reverse procedure-- 
removal of the calling statement and the transfer card and replacement of 
these cards with the transfer card assigned the statement number of the 
deleted subroutine. For example, if equipment subroutine 15 were deleted 
the FORTRAN cards described above would be replaced by the foliowing 
FORTRAN statement: 
Column No.: 4 7 
ISF COTO 7 

One additional change is required in addition to those decribed 
above when an equipment subroutine is incorporated into or removed from 
PACER. This change is necessary due to a feature of PACER which estab- 
lishes the correlation between equipment subroutine names and numbers 
and permits the multiple calling of equipment subroutines for different 
equipments, That is equipments number 5 and 10 may both be pumps 
which can be calculated by subroutine PUMP5. This feature referred to as 
the LPHA listing stores the last three characters of the equipment sub- 
routine name in the position corresponding to its equipment subroutine 
number. The format of the LPHA listing, I2/ (15A3, 10A3), permits the 
inclusion of the first 25 equipments on card 1. If 30 equipments were 
used two cards would be required--25 on the first card and 5 in the first 
15 spaces of the second card. If an equipment subroutine number is not 
occupied, its three characters are left blank. A control card indicating 
the number of equipments for which an LPHA designation should be read 
precedes the LPHA listing and is indicated in the format statement as an 


I2 format. 
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If equipment subroutine COOLR8 were assigned subroutine number 
15, the following entry would be made in column 43-45 of the first card 
of the LPHA listing: 


Column No. 43 - 45 
LRS 





If equipment subroutine PUMPS were designated equipment number 28 the 
following entry would be made in columns 10-12 of the second card of the 
LPHA listing: | 


ОПЕ 10-12 
If, on the other hand, equipment subroutine PUMPS were removed from the 


PACER structure, columns 10-12 of the second card would be left blank. 


Data Deck. The data deck is the means by which the user conveys 
to the PACER executive routine is a detailed discription of the system to 
be simulated. Only one simulation run may be processed for each deta 
deck supplied. However, several simulation runs may be processed 


consecutively by providing consecutive data decks for each run. 


In general, a data deck for a simulation run is composed of 15 
groups of data. However, the first group, the LPHA listing, is provided 
only with the first data deck of a multi-deck execution. All decks after 
the first deck consist of only 14 data groups. The data groups and their 
required order of assembly are listed and discussed below. In addition, 
data decks for the six simulations run to test the feasibility of the moditied 
PACER program are provided as illustrative examples at the end of this 


appendix. (See Figures A3. - A8.). 


LPHA Listing. Format: 12/(15A3,10A3). This data group indicates 
the allowable names of PACER equipment subroutines and provides the 
requisite correspondence between the equipment subroutine number, the 
equipment subroutine name and the equipment number. The first card 
contains in 12 format the number of equipment names to be read from the 
LPHA listing. The remaining cards contain the LPHA listing of the last 
three characters of the allowable equiprnent subroutines. Na redundancy 
is permitted in the three-character names although individual character 
redundancy is permitted. The absolute location of an equipment subroutine 


in the LPHA listing must coincide with the number assigned to the equipment 
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subroutine. That is, the equipment subroutine name placed in the 25th 
location in the LPHA listing must be the last three characters of equipment 


number 25, 


litle Card. Format: 2(12A3). A single card provided for the assign- 


ment of a unique identifying title to each run. 


Control Parameters. Format: 10175. A list of 8 parameters which 
control the processing of PACER simulations. The parameters, in their 


order of listing on the input record, are: 


1. KRUN -~ Run number. 

2. KSETS - Controls intermediate printing of PACER status. 
Following is a listing of allowable values and options each 
provides. 


Value Option 

0 Prints initial conditions of matrices, inter- 
mediate conditions of matrices every KPRINT 
loops and final results of run. 

2 In addition to above, prints values in STRMI 
and STRMO after each caiculation and number 
of stream variable deviations from convercence 
criterias | 

3 Jm addition to above, prints flag lists for 
streams and equipments. 

3. KCLEAN - If set equal to zero, it ceuses PACER to zero 
all matrices and vectors before reading in data. 

4. LOOPS - Maximum number of loops to be permitted in 
iterative calculations. 

5. NOGO - Controls simulation options if convergence is 
not obtained in number of LOOPS permitted. If NOGO = 0, 
the answers obtained during the last loop will be taken as 
correct and the calculation will continue. If NOGO # 0, 
additional trial and error testing of assumed input streams 
will be made until either convergence is obtained, execution 
time expires or calculation is determined to be impossible. 
6. KARDS - If greater than zero, provides a deck of cards 
punched with the stream variables matrix. 

7. KPRINT - Controls number of loops between printout of 
stream variables matrix. 

8. KONV - Not used. Leave blank. 


Dimension Parameters. Format: 1015. This data group contains 
nine parameters which control the allowable dimensions of a run. Any 
change in these variables reauires an identical change to the parameter 
where it appears in the dimension deck. The following is a listing of the 
nine parameters in the order of their recording on the input record: 


1. N2MAX - maximum number of rows permitted in the process 
matrix. 
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2.  N3MAX - maximum number of columns permitted in process 
Matrix. 

3. NEMAX - maximum number of equipments permitted in 
system to be simulated. 

4. NSMAX - maximum number of streams permitted in system 
to be simulated. 

о. NELMAX - maximum number of columns permitted in EN 
Matrix: 

6. NECLMX - maximum number of columns permitted in equip- 
ment control matrix. 

7. NSLMAX - maximum number of columns permitted in SN 
matrix. | 
8. NSCLMX - maximum number of columns permitted in stream 
variables control matrix. 

9. NOCOMP - maximum number of components in each stream. 
Not used other than for information purposes. 

Breoeess Matrix (КРМ), Format: 15, IX2A3, 1X26) 1 Solos 
process matrix contains the topology or configuration of the system to be 
simulated. It is preceded by a single card (Format 15) which indicates 
the number of rows to be read in the process matrix. Each row corresponds 
to an equipment in the system configuration. Equipment numbers need not 
be consecutive but must be in ascending numerical order. The first item 
to be read from the input record (Format I5) is the equipment number. The 
second and third items are the first three and second three characters, 
respectively, of the 6 character equipment subroutine used to calculate 
that equipment. The next 9 items are the input and output streams for that 
equipment. Input streams are listed first and output streams, preceded by 
a negative sign, follow. A word of caution is necessary at this point. One 
must remember that equipment subroutines are written in anticipation of or 
in conjunction with the order of input and output streams. Accordingly, 
any changes in the order of input and output streams for an equipment must 
be checked against the equipment subroutine to insure that revised input/ 
output stream listing is consistent with equipment subroutine. Ifa dis- 
parity exists either the equipment subroutine or the order of input/output 


streams must be revised to bring them into conformity with each other. 


Equipment Parameter Matrix (EN). Format: 5F15.5. The EN matrix 
is preceded by a single control card (Format 15) which indicates the number 
of equipment parameters vectors to 5e read in the EN matrix. The first item 
read in an equipment parameters vector is the equipment number. The 
second item on the input record is not used and is left blank. The remaining 


thirteen items are the equipment parameters for the equipment in question. 





The equipment parameters are unique for each equipment and normally 
represent such parameters as power, heat transfer coefficient, etc. The 
equipment parameters for each equipment employed in the steam power 


system are listed in Figure A9. 


Equipment Control Parameters Matrix (ENC), Format: 5F15.5. The 
ENC matrix is preceded by a single card (Format I5) which indicates the 





number of rows in the ENC matrix. Representative data to be used in the 
ENC matrix wouid normally include allowable differentials across equip- 
ments, costs, etc. Control vectors are only included for those equipments 
which employ them. If no control vectors are employed, the control саға 
should be punched zero or left blank and input/output should proceed to 


next data seament. 


stream Variables Matrix (SN). Format: 5FI5.5,. The SN metrix 15 


preceded by a single card which indicates number of stream vectors to be 
read into program. If contro] variable is zero or blank, the input/output 
will proceed to next data group. The stream variable matrix contains such 
variables as pressure, temperature, enthalpy, total steam flow, etc. 

This location labeling is consistent for each subroutine. The SN vector 


for the steam power system is structured as follows. 





ӨТЕМ, FLAG PRESSURE TEMPERATURE ENTHALPY 
(PSIA) (ШЕ) (BTU/LB) 
SPEC. WT. VELOC, "TOWAL PROV H,O FLOW ОТЕА МРЕ ТО МИ 


(ТВИРТЗ) (FPS) (LB/HR) (LB/HR) (LB/HR) 














INPURITIES FLOW JAIR FLOW | FUEL FLOW |PEUE GAS PEO SV C C И ВИ 
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Stream ControlWMatrix (SNC). Format: 5FI5.5. A*Simgle card 
precedes SNC and indicates number of stream control variable vectors to 
be read. If this value is zero, PACER proceeds to next data group in data 
deck. The SNC matrix contains control parameters such as viscosity, 
volume, etc. The distinction between the variables specified for this 
matrix and that provided for SN are arbitrary and lack discrimination. 
Accordingly, SNC need not be used unless required. It is not employed 


in the steam power cycle. 


Stream Variables Test Vector (DELS). Format: 5FI5.5. This group 


is preceded by a single card (Format I5) which indicates test vectors is 
to be read from DELS matrix. If control card - 0, по test vector is to be 
read and PACER skips to next data group. If contro! card value #0 then 
PACER will read one test vector NSLMAX long. This same test vector 
will be used as test vector for determining convergence of each stream 


variable for all streams. 


Stream Control Variables Test Vector (DELSC). Format: 5FIS.5. 
Preceded by single control card which indicates if test vector is employed. 
If control card value = 0 (Format I5) then no control test vector is employed 
and PACER will proceed to next data group. If control card value #0, 
then test vector is employed and is NSCLMX long. Only one test vector 
is employed for stream control variables for all streams employing control 
variables. This test vector will check convergence of controi variables 


for all streams employing control variables. 


Preferred Stream Numbers (KPS). Format: 1015. KPS is preceded by 


a control card (Format I5) which indicates number of preferred streams. 

If control card value equals zero, then PACER proceeds to next data group. 
If control card 40, PACER reads numbers of preferred streams, i. e. those 
streams which will be assumed known first in a trial and error calculation 
to determine order of calculations. Normally these should be recycle or 
other internal streams since feed and product streams are assumed known 


in trial and error calculations. 


Extra Equipment Vectors (AEN). Format: 5FI5.5. AEN is preceded 
by a control card. If control card value (Format I5) equals zero, PACER 
will proceed to next data grouo. If contro! value 40, PACER reads into 


AEN extra equipment parameters for number of equipments indicated on 
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control card. This matrix is only used if insufficient space is available 


in EN and ENC. Length of vector is arbitrarily established at 40. 


Order of Equipment Calculation (NELIST). Format: 1015. NELIST 


is preceded by a control card which indicates number of equipments con- 
tained in NELIST matrix. If the order of calculation is not specified, 

the control card value (Format 1015) equals zero and PACER proceeds to 
next data group. Each card consists of five pairs of data. The first 
value of a pair represents the equipment number. The second value 
represents the calculational mode for that equipment. If mode indicator 
equals one, the equipment will be calculated directly. If the mode 
indicator equals two, the equipment will be calculated as an equipment 
in an iterative list of equipment, Finally, if mode indicator equals three, 
the equipment will be calculated as last item in list of equipment to be 


calculated in an iterative mode, 


Selected Streams Vector (KSPRNT). Format: 1615. Contains as the 
first item the number of stream variable vectors to be printed each itera- 
tion. If it is not desired to print any stream variable vectors each itera- 
tion the card should be left blank. The stream numbers for which the 
stream variables are to be printed are listed after the first item. A maxi- 


mum of 15 stream variable vectors may be printed each iteration. 
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Appendix B 


RESULTS OF TESTS 


Format and Interpretation of Results. As discussed in Chapter V, 
six simulation tests were developed to test the validity of the revisions 
made to the original PACER program, the efficacy of PACER as a simulation 
device and the potential usefulness of PACER as an instructional device 
in a chemical engineering curriculum. The results of those simulation 


tests are presented as Figures Bl. - B6. 


The final values of the stream variables matrix and the equipment 
parameters matrix are presented for each run. In addition for runs 01 and 
03 extensive diagnostics of the PACER solution process are presented 
for analysis by the user. This is a PACER option. The stream variables 
are consistent for all streams in the stream variables matrix. That is, 
stream variable 5 will represent enthalpy for allstreams. However, the 


enthalpy for all streams will not necessarily be the same. 


For the steam power system developed for this project, the following 
stream variable arrangement was decided on and should be used as a 
guide in interpreting values of stream variables in the stream variables 
matrix. Note that there are five columns in each stream. This dimension 
may be parametically varied by the user by adjusting dimensions of SN 


matrix. 


STRM. FLAG PRESSURE | TEMPERATURE ENTHALPY 
(PSIA) (Әс) (BTU/LB) 
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